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TECHNICAL MEMORANDUM 


STS-35 SCRUB 3 HYDROGEN LEAK ANALYSIS 
INTRODUCTION 


Although considerable analysis had been performed on the previous STS-35 scrub and 
tanking test data, the wide range of operating conditions and valve configurations tested during 
scrub 3 was not previously available. With the formation of the STS-35 leak investigation team, 
an intensive analysis effort was undertaken to characterize the leak and determine its location. 
This report describes the analysis in roughly the chronological order in which it was performed. 

The results of this analysis effort were provided to the leak investigation team as they 
became available. The results presented here sometimes differ slightly from those provided to 
the team during the investigation due to error corrections and refinements. However, the conclu- 
sions drawn from the results remain unchanged. 

During the STS-35 scrub 3 leak isolation test, the behavior of the aft compartment H 2 
concentration seemed to indicate a leak in the engine 2 prevalve (PV5). The decrease in H 2 
concentration during drain associated with opening the engine 2 prevalve was particularly indit- 
ing. However, since the number of different leaks was unknown, all possible locations were con- 
sidered in the analysis. 


SUMMARY 


Some of the more important conclusions drawn from this analysis were: 

1. The leak did not exist at ambient temperature. 

2. The engine 2 prevalve was the most likely leak location. 

3. At least 80 percent of leakage came from the engine 2 prevalve. 

4. The scrub 2 leak area was twice that of scrub 3 and consistent with the known engine 
3 prevalve detent cover seal leakage. 

5. Leak area changes cannot be inferred from concentration changes without employing an 
analysis similar to that used here. 

MODEL METHODOLOGY 

The basic approach taken was to calculate an aft compartment H 2 concentration based on 
the measured pressures, temperatures, and valve configurations and evaluate this prediction by 
direct comparison with the mass spectrometer data. 


Prior to STS-35 scrub 3, H 2 flow rates had been calculated for an adiabatic real gas 
expansion over a range of initial pressures and temperatures. These data were converted to a 
table which provided relative flow rate as a function of the leak site pressure and degrees of 
subcooling. 

Aft compartment H 2 concentration was calculated by assuming that the H 2 was uniformly 
mixed in a single volume and using the transient mass conservation equation: 

d Mh 2 


where. 


Y c A*N 2 

M N2 = mass of N 2 in aft compartment 


mn, = H 2 leak rate 

*in 


" IH 2ou. = Yc WN 2 

wn 2 = aft compartment purge mass flow rate 
Y c = aft compartment H 2 concentration. 


While this approach preserves the basic time constant of the aft compartment purge flow 
(about 90 s), the calculated H 2 concentration begins to respond immediately to a change in leak 
flow rate. The discrete H 2 transport delay from the leak site to the vent door is smeared over the 
time constant. In addition, the predictions do not include time delays associated with the mass 
spectrometer H 2 concentration measurement. 


SCRUB 3 MEASURED DATA 


Rather than using the measured manifold pressure, a value calculated from measured 
ullage pressure and estimated liquid height was used. Both the measured and calculated manifold 
pressures are shown in figure 1. 

The measured manifold temperature is compared with the saturation temperature cor- 
responding to the calculated manifold pressure in figure 2. In this study, all saturation tempera- 
tures were determined from the pressure by the approximation: 

TSAT(DEGR) = 9.96915 (psia)0.2«8438 +15.0744 . 

For leak rate predictions, temperature data were typically corrected for offsets by compar- 
ing the measurement to calculated saturation temperatures at appropriate times. 
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Engine inlet temperatures are compared with the calculated manifold saturation tempera- 
ture in figures 3, 4, and 5. Engine inlet pressures are compared with the calculated manifold pres- 
sure in figures 6, 7, and 8. 

The midrange 5,000 parts per million (p/m) aft compartment H 2 concentration measure- 
ments, GGDR2510T, was used to compare with all predictions. However, as shown in figure 9, 
there were occasional peaks above 5,000 p/m which were clipped by this measurement. 

An event timeline was constructed which provides vehicle and facility valve actuation 
times and is shown in table 1 . All data plots were made as seconds from the Space Transporta- 
tion System (STS) data base Greenwich mean time (GMT) reference time of 90:260:21:47:00.0. 


INITIAL MANIFOLD/ENGINE LEAK PREDICTION RESULTS 


Initially, leak rate predictions were performed for the manifold based on the calculated 
manifold pressure and measured manifold temperature (V41T1428A). Predictions for each engine 
based on measured engine inlet pressure and measured inlet temperature were also performed. 
These results are compared with the 5,000 p/m range measured data in figures 10 to 13. These 
predictions did not include any heat leak effects on the leak site temperature. 

The comparisons show the manifold prediction correlates much of the data except for the 
drain time period after 25,000 s. The engine 1 and engine 3 predictions show poor correlation. 
However, the engine 2 comparison shows good correlation during the drain period. 


SPECIAL TEST PREDICTION RESULTS 


After the initial manifold and engine predictions were performed, attention was focused on 
four specific time periods when unusual conditions and H 2 concentration behavior was observed. 


A. Recirculation Pump Backspin 

The first of these special conditions considered occurred between 9,000 and 10,000 s 
when, with engines 2 and 3 prevalves closed, the recirculation (recirc) valves were opened and 
engine 1 recirc began. The aft compartment H 2 p/m dropped as shown in figure 14. This sequence 
caused the engine 2 and engine 3 recirc pumps to backflow as shown by the recirc pump speed 
data in figure 15. 

The engine 1 inlet pressure and temperature are shown in figures 16 and 17. Engines 2 
and 3 inlet temperatures remained offscale high. The manifold temperature also increased during 
this time. 

The recirc pump backspin was clearly the result of vapor flow from the warm engine inlet 
lines. Based on the recirc pump acceleration and estimated pump delta p, a vapor flow rate of 0.02 
lb/s per pump was calculated. 
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Predicted leak rates were calculated for a leak site temperature corresponding to sat- 
urated liquid and saturated vapor when the recirc valve was open and the manifold temperature 
when closed. The comparisons in figure 18 show good agreement between the measured concen- 
tration and the saturated vapor prediction. However, a 120-s time delay was found which could 
not be totally explained. 

The conclusion drawn from this analysis was that the vapor flowing out of recirc pumps 2 
and 3 inlet produced vapor at the leak site. Since this occurred during reduced fast-fill and the 
flow patterns in the manifold were unknown, the leak could have been in the manifold, the 
prevalves, or the 17-in line. 

B. Prevalve Opening at 14,000 s (GMT 261:01:40) 

After engine 3 recirc was terminated and the recirc valves were closed, all three prevalves 
were opened and the p/m level dropped. Subsequent closing of the prevalves resulted in the p/m 
data returning to their previous level. During this period, the engine 1 bleed valve was closed and 
the inlet line dry, the engine 2 bleed valve was closed and the inlet line saturated, and engine 3 
remained subcooled. The manifold temperature increased when engines 1 and 2 prevalves were 
opened. 

Both engines 1 and 2 were probably producing vapor and a recirculating flow in their inlet 
lines. As in the previous special test analysis, a prediction was made using the manifold tem- 
perature when the engine 2 prevalve was closed and saturated liquid or vapor when the engine 2 
prevalve was open. The comparison between the saturated vapor prediction and measured con- 
centration shown in figure 19 indicated good agreement with about a 30-s time delay. 

The conclusion drawn from this analysis was that the leak was probably not in the engine 
3 prevalve, but further discrimination of the low pressure components was not possible. 


C. Initial Fast-Fill and Transition to Reduced Fast-Fill 

When the manifold flow rate was reduced from the 8,300 gallons per minute (gal/min) 
fast-fill rate to the 930 gal/min reduced fast-fill rate, the measured p/m dropped. About 10 min 
later, the concentration returned to the previous value. 

The primary effect of manifold flow on the leak rate is the effect on leak site temperature. 
As a first approximation, the leak site temperature should be higher than the bulk fluid tempera- 
ture by an amount proportional to the inverse of the flow rate. A prediction was made assuming a 
leak site temperature equal to the manifold temperature plus 0.6° during fast-fill and plus 6° 
during reduced fast-fill. The leak site temperature was arbitrarily returned to the manifold tem- 
perature at 3,400 s. The rise in leak rate at this time was thought to be the result of chill down 
and filling of the volume at the leak site and could not be reasonably modeled. Prior to 1,600 s, 
the predicted leak rate was set to zero. Figure 20 compares this prediction with the measured 
concentration. 

Three conclusions were drawn from this analysis. The predicted 600 p/m from a saturated 
vapor leak was not observed in the data. Even a 200-p/m ambient temperature H 2 leak should 
have been detectable by the mass spectrometer. As shown in figure 21, the measured concentra- 
tion was less than 30 p/m. This low measured concentration relative to the predicted value 


4 



confirmed that the effective leak area existed only when the leak site was chilled to operating 
temperatures and that an ambient temperature leak did not exist. 

The second conclusion was that the manifold temperature and assumed heat leak effect 
was consistent with the measured concentration. 

The third conclusion was that, with the prediction available, an ullage pressure effect on 
the leak rate could be seen in the measured concentration. 


D. Prevalve Opening During Drain 

During drain (starting at about 25,000 s), each engine's recirc pump was stopped and its 
prevalve opened. While there was no effect on the p/m data for engines 1 and 3, when engine 2's 
prevalve was opened the p/m level dropped. The engine 2 effect repeated three times. 

As described in case b, each engine was probably producing vapor in its inlet line when its 
prevalve was open. As previously mentioned, the concentration data behavior was consistent 
with the engine 2 inlet temperature. Since the 5,700 gal/min drain flow produced a flow velocity of 
about 8 ft/s and there was no effect on the manifold temperature measurement and since the 
17-in line is essentially horizontal, this vapor was very unlikely to rise above the manifold. Thus, 
these data excluded the 17-in line and the disconnect as possible leak sites. Further, it also 
excluded the engines 1 and 3 prevalves, since they produced no effect. These data, together with 
the recirc pump backspin data, indicated that the engine 2 prevalve, PV5, was the most likely 
leak location. 


PV5 SCRUB 3 PREDICTION RESULTS 

While the special test analysis had provided many important conclusions and given much 
insight into the leak behavior, it was felt that a unified model containing all the logic developed to 
date was required. This unified model would permit verification of all of the assumptions and logic 
through comparison with all of the concentration data. 

The unified model consisted of a set of simple rules for determining the leak site tempera- 
ture. The leak site temperature was taken to be the manifold temperature when the engine 2 
prevalve was closed, otherwise it was taken to be the engine 2 inlet temperature. A heat leak 
effect was approximated by adding (930 gal/min/manifold flow) to this temperature. This may 
overestimate the heat leak effect but was judged acceptable. The manifold flow used is shown in 
figure 22. 

To describe the effect of recirc pump vapor backflow, special logic was used to override 
the leak site temperature and instead used a saturated vapor leak rate. The logic used was that 
vapor backflow would occur when the recirc valve was open, when the recirc pump was off, when 
the prevalve was closed, and when the engine inlet temperature was superheated. 

When this logic was initially implemented, the model predicted the concentration reduc- 
tion previously analyzed. However, as shown in figure 23, the model also predicted backflow and 
concentration reduction at about 11,000 s which was previously unnoticed. In addition, when ini- 
tially implemented, the logic was applied equally to till three engines. Engine 1 was predicted to 
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have vapor backflow between 21,100 and 12,600 s and between 12,900 and 14,100 s while no 
concentration reduction occurred in the data. Apparently, the local flow field transported the vapor 
from the engine 2 recirc inlet to the leak site but did not transport the vapor from recirc pump 1. 
This result was considered to be additional evidence of an engine 2 prevalve leak since each 
recirc pump is located directly in front of its prevalve. The final version of the model only applies 
the backflow logic to engine 2. 

The model was run on a PC. A listing of the program is provided in the appendix. The 
model and all test data used in this analysis are available from the author on a PC floppy disk. 

The model prediction is compared with the measured H 2 concentration in figure 24 for the 
full duration of the data and in figures 25 to 29 with an expanded time scale. Also shown in these 
figures are event bars which indicate when valves are open or the recirc pumps are on. Figure 30 
shows the predicted leak site subcooling along with the manifold subcooling. 

Two deviations of the prediction from the measured data are shown in figure 26. The first, 
at 11,400 s, occurs when closure of the recirc valves terminates vapor backflow throughout the 
engine 2 recirc pump and the model switches from a saturated vapor leak rate to a leak based on 
the subcooled manifold temperature. This deviation is thought to be caused by the time required 
to condense the accumulated vapor volume in the prevalve which is not modeled. A similar time 
delay was found when vapor backflow was terminated at 9,400 s. 

The second deviation occurs at 12,700 s when the engine 2 prevalve is opened, and the 
model predicts a concentration decrease corresponding to the near saturated engine inlet tem- 
perature. The cause of this deviation is unknown. Similar prevalve openings between 14,000 and 
15,000 s and during drain produced a concentration decrease as predicted. 

To determine if all of the scrub 3 leak could be explained by one PV5 leak, several com- 
binations of a PV5 leak and a manifold leak were predicted. As previously described and shown 
in figure 10, a pure manifold leak prediction does not include any effect of recirc pump vapor back- 
flow or prevalve opening. Thus, combining the PV5 leak with a manifold leak tends to reduce the 
magnitude of the predicted concentration reduction for these effects. Figure 3 1 shows a prediction 
for an 80-percent PV5 leak combined with a 20-percent manifold leak. Further increases in the 
manifold leak fraction tends to raise the minimum concentration during the vapor backflow and 
PV5 open times to a value above the data. Therefore, the 20-percent manifold fraction was 
judged to be the highest reasonable fraction which could exist and most of the leak was in PV5. 


PV5/PV6 SCRUB 2 PREDICTION RESULTS 

To further verify the scrub 3 model, predictions were made for scrub 2. Figure 32 shows 
the calculated manifold pressure, and figure 33 shows the manifold flow rate. Initially, the scrub 3 
model was applied to scrub 2 without any modification. Surprisingly, the predicted concentration 
was almost exactly half the measured data. Although the PV6 detent cover seal was a known 
additional leak source for scrub 2, the measured scrub 2 concentration was only slightly higher, 
and this result was totally unexpected. 

Subsequent scrub 2 predictions modeled the PV6 leak as a separate second leak with the 
same leak area and same logic as developed for PV5. The final scrub 2 prediction is compared 
with the measured concentration in figure 34. The leak site subcooling and manifold subcooling 
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are shown in figure 35. The reasonably good agreement further verified the model but, more 
importantly, provided further evidence that most of the scrub 3 leak was due to the one PV5 leak. 


These results clearly showed that the changes in concentration data between scrub 2 and 
scrub 3 could not be used to infer changes in leak area. The effects of the different valve config- 
urations and different manifold flows could only be included by applying the model analysis. The 
scrub 3 leak area was actually half that of scrub 2. 


PV5/PV6 TANKING TEST 1 PREDICTION RESULTS 


This tanking test was primarily concerned with the orbiter/extemal tank disconnect leak- 
age. High external H 2 concentrations resulting from the disconnect leakage prevented high fill 
rates. As shown in figure 37, the manifold flow rate was low for most of the test. 

The model prediction H 2 concentration is compared with the measured data in figure 38. 
This prediction used a leak area 83 percent of that used for scrub 2. There is considerable uncer- 
tainty in the actual amount of leak area reduction since significant concentrations were obtained 
for only two brief periods. 


PV5/PV6 SCRUB 1 PREDICTION RESULTS 


The scrub 1 manifold flow conditions, shown in figure 41, were very similar to tanking 
test 1. Engine recirculation was not performed. The model prediction H 2 concentration is 
compared with measured data in figure 42 for a leak area of 69 percent of that used for scrub 2. 
Again, considerable uncertainty exists in the amount of this leak area reduction. However, since 
some additional H 2 leak flow into the aft compartment from the disconnect leak should have 
existed during scrub 1 and tanking test 1, there does appear to be a trend of increasing leak area 
with each loading. 


CONCLUSIONS 


The leak isolation procedures used during scrub 1 and tanking test 1 were designed to 
further understand the disconnect leak and did not produce useful data relative to an aft compart- 
ment leak. With sufficient amounts of the right test data, detailed analysis and modeling can 
determine the leak location and behavior. Although the results of the STS-35 analysis were not 
always sufficiently timely to significantly contribute to the leak investigation team s activities, 
future leak investigation analyses can build on the experience and analysis tools developed here. 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline. 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline (continued) 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline (continued) 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline (continued) 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline (continued) 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline (continued) 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline (continued) 
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Table 1. STS-35 scrub 3 LH 2 leak isolation timeline (continued) 
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Figure 1. STS-35 S3 manifold pressure. 
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Figure 2. STS-35 S3 manifold temperature 
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Figure 3. STS-35 S3 engine 1 inlet temperature. 
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Figure 4. STS-35 S3 engine 2 inlet temperature. 
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Figure 7. STS-35 S3 engine 2 inlet pressure. 
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Figure 14. STS-35 S3 aft compartment H 2 concentration. 








CD LO CD ID CD 

K VO VO ID ID 


y ist 


32 


45 . 


o 

n 





o 

o 

o 

° cn 

• 

• 

• 

o 

Lf) 

o 


CO 

CO 


Figure 16. STS-35 S3 engine 1 inlet pressure. 
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Figure 17. STS-35 S3 engine 1 inlet temperature. 
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Figure 18. STS-35 S3 predicted aft compartment H 2 concentration. 
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Figure 19. STS-35 S3 predicted aft compartment H 2 concentration. 












Figure 20. STS-35 S3 fast-fill/reduced-flow transition. 
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Figure 21. STS-35 S3 low range aft compartment H 2 concentration. 
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Figure 22. STS-35 S3 manifold flow. 
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Figure 23. STS-35 S3 engine 2 recirculation pump vapor. 
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Figure 27. STS-35 S3 predicted PV5 leak. 
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Figure 29. STS-35 S3 predicted PV5 leak. 
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Figure 30. STS-35 S3 leak site subcool temperature. 







Figure 31. STS-35 S3 80-percent PV5 + 20-percent 17-in line leak. 









Figure 32. STS-35 S2 manifold pressure. 
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Figure 34. STS-35 S2 predicted versus actual aft H 2 concentration. 
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Figure 35. STS-35 S2 leak site subcool versus manifold subcool. 
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Figure 36. STS-35 T1 manifold pressure. 
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Figure 37. STS-35 T1 flow. 
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Figure 39. STS-35 T1 leak site subcool versus manifold subcool. 




Figure 40. STS-35 Si manifold 
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Figure 41. STS-35 SI flow. 



ACTUAL PREDICTED 



Figure 42. STS-35 SI predicted versus actual aft H 2 concentration. 




LEAKSITE SUBCOOL MANIFOLD SUBCOOL 



Figure 43. STS-35 SI leak site subcool versus manifold subcool. 







APPENDIX 


PRECEDING PAGE BLANK NOT FILMED 



c Routine to calculate Aft Comp PPM based on subcool inq and oressure 
c STS35S3 PV5 leak model 
c 


c 


c 


c 

c 

c 


c 


parameter ( 

vol 

= 

4550 

• ) 

I aft comp volume (ft3) 

parameter ( 

volin3= 

Vol*1728 . ) 

! " " (in3) 

parameter ( 

purge 

= 

235. 

) 

! aft comp purge (lb/min) 

parameter ( 

w 

= 

1728 

.♦purge/. 072 ) 

! purge flow (scim) 

parameter ( 

nrecl 

= 

2 ) 

! recircl 


parameter ( 

nrec2 

= 

3 ) 

! recirc2 


parameter ( 

nrec3 

= 

4 ) 

! recirc3 


parameter ( 

npvl4 

= 

5 ) 

!pvl4 


parameter ( 

npvl5 

= 

6 ) 

! pvl5 


parameter ( 

npvl6 

= 

7 ) 

!pvl6 


parameter ( 

npv4 

= 

8 ) 

!pv4 


parameter ( 

npv5 

— 

9 ) 

Ipv5 


parameter ( 

npv6 

= 

10) 

!pv6 


parameter ( 

npv22 

= 

ID 

! pv22 


parameter ( 

npvl3 

= 

12) 

! pvl3 


parameter ( 

npvl2 

= 

13) 

! pvl2 


parameter ( 

nbldl 

= 

14) 

! elbld 


parameter ( 

nbld2 

= 

15) 

! e2bld 


parameter ( 

nbld3 

= 

16) 

! e3bld 


real al(34) 






real a2(22) 






real event (17) 





real tempin(3) 





LEAK RATE TABLE 





data al/32. 

9 0 • / “ 

■15 

.35, 

2.16, -11.35, 

2.12, -7.35, 1.82, 


-5.35, 1.62, -4.35, 1.5, -3.35, 1.37, -2.35, 1.24, 
-1.35, 1.14, -.35, 1.03, 0.0, 1.0, .1, .28, .5, .274, 

.65, .27, 4.65, .25, 14.65, .22, 24.65, .21/ 

data a2/20 . 0 , 0., 20., .35, 25., .50, 30., .60, 

35., .72, 40., .81, 45., .91, 50., 1.0, 

55., 1.09, 60., 1.18, 70., 1.28 / 


TSAT(psia) = 9. 969152*(psia)**. 2843842 + 15.07444 
c 
c 

open(unit=l, file=' d:\sts35s3\mandt.bin* , status^old* , 

* f orm= 1 unformatted 1 , access= 1 direct 1 ) 
c 

open (unit=2 , f ile= 1 d : \sts35s3\manpr . cal 1 , status= • old 1 , 

* f orm= 1 unformatted 1 , access= 1 direct 1 ) 
c 

open (unit=3 , f ile= 1 d : \sts3 5s3\manpred . bin 1 , f orm= 1 unformatted 1 , 

* status= 1 unknown 1 , access= 1 direct 1 , recl=4 0 ) 
c 

open (unit=7 , f ile= 'd:\sts35s3\t\events.evn 1 , status= 1 old 1 ) 
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c 


★ 

* 

* 



status= , old l , 
status= , old' , 
, status^ 1 old 1 


/ 


c 


c 


c 

c 

c 


c 

c 

c 

200 


c 

c 

c 


write (5, • (a) ' ) ' Enter Tsat Tolerance ' 

read (5,*) tolinp 

read ( 1 ) timb,tempb 
read (2 ) tim2,prsl 
read (8 ) timel, temel 
read (9 ) time2,teme2 
read (10) tiroe3,teme3 


tim = tim2 
timxv = tim 
tprint= tim 
prs2 = prsl 
ppml = 0 . 
ppm2 = 0 . 
ppm3 = 0. 

MAIN TIME LOOP 


continue 
tim=tim+l . 
tol = tolinp 

SET MANIFLOD FLOW RATE 


! Time step fixed at 1 sec 


if ( tim .It. 2590 ) then 
flow = 8300. 

else if( tim .It. 17004. ) then 
flow = 930. 

else if ( tim .It. 18392. ) then 
flow = 600. 

else if ( tim .It. 20406. ) then 
flow = 930. 

else if( tim .It. 25801. ) then 
flow = 600. 

else if ( tim .It. 35000. ) then 
flow = 5700. 

end if 


205 continue 

if ( tim .It. timb ) goto 210 
tima=timb 
tempa=tempb 

read(l,end=900) timb,tempb 
tprint=tim-l . 
goto 205 


Update manifold temp 
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c 


210 

1 


c 

215 
c 

c 

c 

c 

c 

c 

c 

216 
c 

c 1 
c 

315 


316 


415 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


416 


218 


c 

c 


219 

220 


continue 

if( timb-tima .gt. .01 ) 

tempm=tempa+ (tim-tima) * (tempb-tempa)/ (timb-tima) 
continue 

if( tim .It. timel ) goto 216 
timelp=timel 
temelp=temel 

read (8 , end=900) timel, temel ! Update engine 1 temp 

tprint=tim-l. 
goto 215 

continue 

if( timel-timelp .gt. .01 ) 

tempin (1) =temelp+ (tim-timelp) *(temel-temelp) / (timel-timelp) 
continue 

if( tim .It. time2 ) goto 316 
time2p=time2 
teme2p=teme2 

read (9 , end=900) time2,teme2 ! Update engine 2 temp 

tprint=tim-l. 
goto 315 

continue 

if( time2-time2p .gt. .01 ) 

tempin (2 ) =teme2p+ (tim-time 2 p) *(teme 2 -teme 2 p)/(time 2 -time 2 p) 

continue 

if( tim .It. time3 ) goto 416 
time3p=time3 
teme3p=teme3 

read ( 10 , end=900) time3,teme3 I Update engine 3 temp 

tprint=tim-l. 
goto 415 

continue 

if( time3-time3p .gt. .01 ) 

tempin ( 3 ) =teme3p+ (tim-time3p) * (teme3-teme3p) / (time3-time3p) 
continue 

if( tim .It. timxv ) goto 219 

read(7, fmt=*,end=900) timxv, event ! Update events 

read ( 7 , f mt=* , end=9 0 0 ) timxv ! Look ahead for next time 

backspace (7) 
tprint=tim-l. 
goto 218 

continue 

if (tim. IT. tim2) goto 230 

timl=tim2 

prsl=prs2 
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c 


r ead ( 2 , end=9 0 0 ) tim2,prs2 
tprint=tim-l . 
goto 220 


i Update pressure 


230 


c 

c 

c 

c 


c 

c 

c 


1 

2 

3 


c 

c 

c 

c 


c 

c 


c 

c 

c 


c 

c 


continue 

if( tim2-timl .gt. .01 ) 

prs=prsl+ (tim-timl) * (prs2-prsl)/ (tim2-timl) 


CHECK FOR VAPOR AT PV5 INLET FROM ENG 2 RECIRC PUMP INLET 
vapor = 0.0 

if ( event (npvl5) .gt. .1 .and. 

event (nrec2) .It. .1 .and. 
event (npv5 ) .It. .1 .and. 

tempin(2) .gt. 1. ) vapor = vapor+1. 

SET LEAK TEMPERATURE 

if( event ( npv5 ). gt . .1 ) then 

temp = tempin(2) 
else 

if( vapor .gt. .1 ) then 

temp = .1 ! Assume sat vapor leak 

else 

temp = tempm 
end if 
end if 

temp * temp+1. *930 ./flow 
GET LEAK FLOW 

call intrp(al,temp-tol, factla) 
call intrp(al, temp, fact lb) 
call intrp(al, temp+tol, factlc) 
call intrp(al, tempm, fact lm) 

call intrp(a2 ,prs, fact2) 

wleakl = factla*fact2 

wleak2 = factlb*fact2 

80% PV5 , 20% manifold case 

wleak2 = . 8 *factlb*fact 2 +. 2*factlm*fact2 

wleak3 = factlc*fact2 

if (tim.lt. 2122.) then ! Ambient temp leak 

wleak3 = .21*fact2 

wleak3=wleak3*sqrt( (TSAT(prs)+temp+tol)/540. ) 

wleak2=wleak3 
wleakl=wleak3 
end if 

GET PPM 
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c 

ppml = ppml+(wleakl-ppml*l.e-6*w)*l.e6/volin3/60. 
ppm2 = ppm2+(wleak2-ppm2*l.e-6*w) *l.e6/volin3/60. 
ppm3 = ppm3+(wleak3-ppm3*l.e-6*w) *l.e6/volin3/60. 
c 

if( tim.gt.O. .and. tixn .ge. tprint ) then 
tprint = tim+10. 
write (3) tim,ppml,ppm2 ,ppm3 , 

1 wleakl,wleak2,wleak3, 

2 temp, flow, vapor 
end if 

c 

goto 200 
c 

900 close (unit=l) 
close (unit=2) 
close (unit=3) 
close (unit=4) 
c 

end 

c 

c 

SUBROUTINE INTRP(A,X,Y) 

C 

C General purpose interpolation routine. A() is the array to 
C interpolate, X in the independent variable, and Y is the 
C returned dependent variable. The A() array should be 
C configured as follows: 

C 

C a(l) Real - number of x and y entries (n = 2 * x/y pairs) 
c a (2) Real - zero 

C a (3) Real - XI 

C a (4) Real - Y1 

C a (5) Real - X2 

C 
C 

C a(n+2) Real - Ym where m=n/2 
C 

C Note that both a(l) and a (2) are reset by INTRP after the first 
C call, and should not be changed by the calling program 
C 

REAL A ( * ) ,A1,A2 
INTEGERM 11,12 
EQUIVALENCE (A1 , II) , (A2 , 12 ) 

C 

A1=A(1) 

A2=A (2 ) 

IF (12 .NE. 0) GOTO 20 ! check if initialized 

I1=A1+1 ! reset to integer 

A(1)=A1 ! and store 

12=3 ! set current index at bottom 

C 

20 IF (X-A (12 ) ) 30,40,50 
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30 


40 

45 


C 


50 


C 


60 


IF ( 12 . EQ. 3 ) GOTO 40 
12 = 12-2 
GOTO 20 


1 at bottom, use first value 
! else backup 
! and try again 


Y=A (12+1) 
A ( 2 ) =A2 
RETURN 


! use current value 
! restore current index 
! and done 


IF (I2.EQ.il) GOTO 40 J at top, use last value 
D1=A (12+2 ) -X 

IF(Dl.LE.O) GOTO 60 ! not far enough 

Y=A (12+3 ) - Dl* (A(I2+3) -A (12+1) )/ (A (12+2) -A (12) ) 

GOTO 45 


12 = 12+2 
GOTO 20 
END 


move foward 
and try again 
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